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1 Problem 

Some biological systems consist of two "phases" of nearly square fiber bundles of differing 
thermal and electrical conductivities. Consider a circular region of radius a near a corner of 
such a system as shown below. 




Phase 1, with electrical conductivity o\, occupies the "bowtie" region of angle ±a, while 
phase 2, with conductivity o" 2 <C 0"i, occupies the remaining region. 

Deduce the approximate form of lines of current density j when a background electric 
field is applied along the symmetry axis of phase 1. What is the effective conductivity o of 
the system, defined by the relation / = aA<f) between the total current / and the potential 
difference A<p across the system? 

It suffices to consider the case that the boundary arc (r = a, \9\ < a) is held at electric 
potential 0=1, while the arc (r = a, 7r — a < \9\ < it) is held at electric potential <fi = — 1, 
and no current flows across the remainder of the boundary. 

Hint: When <j 2 <C a±, the electric potential is well described by the leading term of a 
series expansion. 

2 Solution 

The series expansion approach is unsuccessful in treating the full problem of a "checkerboard" 
array of two phases if those phases meet in sharp corners as shown above. However, an 
analytic form for the electric potential of a two-phase (and also a four-phase) checkerboard 
can be obtained using conformal mapping of certain elliptic functions pQ. If the regions of 
one phase are completely surrounded by the other phase, rather lengthy series expansions 



for the potential can be given |2j. The present problem is based on work by Grimvall [3] and 
Keller 0|. 

In the steady state, the electric field obeys V x E = 0, so that E can be deduced from a 
scalar potential <p via E = — V0. The steady current density obeys V ■ j = 0, and is related 
to the electric field by Ohm's law, j = crE. Hence, within regions of uniform conductivity, 
V ■ E = and V 2 </> = 0. Thus, we seek solutions to Laplace's equations in the four regions of 
uniform conductivity, subject to the stated boundary conditions at the outer radius, as well 
as the matching conditions that <fr, E\\, and j± are continuous at the boundaries between the 
regions. 

We analyze this two-dimensional problem in a cylindrical coordinate system (r, 9) with 
origin at the corner between the phases and 8 = along the radius vector that bisects the 
region whose potential is unity at r = a. The four regions of uniform conductivity are labeled 
/, 77, /// and IV as shown below. 



dtydr = 



<\>=-1 




d$/dr = 

Since jj_ = j r = oE r = —ad<f)/dr at the outer boundary, the boundary conditions at 



r = a can be written 

Mr = a) = 1, (1) 
d</> n (r = a) d(p IV {r = a) = Q 

dr dr 

<frin(r = a) = -1. (3) 

Likewise, the condition that j± = jg = aE e = —(a/r)d<f)/d9 is continuous at the boundaries 
between the regions can be written 

d(J) I (6 = a) d<t> n (0 = a) 

ai — m — = ° 2 — w — ' (4) 

d(f>m{9 = 7r - a) d(j} H {9 = ir — a) , . 

ai 86 = ° 2 M ' (5) 

etc. 

From the symmetry of the problem we see that 

(6) 
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0(^-0) = -0(0), (7) 

and in particular 0(r = 0) = = 0(0 = ±ir/2). 

We recall that two-dimensional solutions to Laplace's equations in cylindrical coordinates 
involve sums of products of r ±k and e ±lk0 , where k is the separation constant that in general 
can take on a sequence of values. Since the potential is zero at the origin, the radial function 
is only r k . The symmetry condition © suggests that the angular functions for region / be 
written as cos k6, while the symmetry condition (J7J) suggests that we use sin/c(7r/2 — |0|) in 
regions II and IV and cos k{n — 0) in region III. That is, we consider the series expansions 

7 = ]TA fe r fe cos £0, (8) 
B k r k sink {6^ , (9) 

-J2A k r k cosk(n-6). (10) 

The potential must be continuous at the boundaries between the regions, which requires 

Ak cos ka = Bk sin k ^— — c^J . (11) 

The normal component of the current density is also continuous across these boundaries, so 
eq. (JD) tells us that 

<j\Ak sin ka = oiBk cos k ^— — a^j . (12) 
On dividing eq. (fT2^) by eq. (fTTj) we find that 

tan ka = — cot k ( — — a ) . (13) 



<j>in 



a i V2 

There is an infinite set of solutions to this transcendental equation. When o^/^i ^ 1 wc 
expect that only the first term in the expansions JBJ-flHJ) will be important, and in this case 
we expect that both ka and k(ir/2 — a) are small. Then eq. (J 13)) can be approximated as 

ka « a2/(Jl (U) 

and hence 

e M <<: L (15) 

Equation (jllj) also tells us that for small ka, 



7T 



A k « ^- - a J . (16) 

Since we now approximate 0/ by the single term A k r k cos k9 ~ A^r : , the boundary 
condition JT} at r = a implies that 

A k « 3 , (17) 
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and eq. (|TB|) then gives 

The boundary condition (J2J) now becomes 



= kB k a sin « — i (19) 

V2 / a(f-a) v ; 

which is approximately satisfied for small k. 

So we accept the first terms of eqs. (JH|)- (|lUj) as our solution, with k, A k and B k given by 
eqs. (USD, (JO and (EHJ. 

In region J the electric field is given by 

E r = — ^ « -k^^coskO « -fc— 7T- , (20) 
1 (90/ , r' " 1 

Thus, in region /, E e /E r ?s k9 <C 1, so the electric field, and the current density, is nearly 
radial. In region // the electric field is given by 

Er = ^»- kT ^ sink (l-e)*- k ^l^±, (22) 



E e = — ^-^k^smk9^k 2 9-^-. (21) 



dr ka k (^-a) \2 J a k § - a 



2 "V v ^2 



* = (23) 



r 89 ka k (\-a) \2 J a k {\ - a] 

Thus, in region II, E r /Eg « k(n/2 — 9) <C 1, so the electric field, and the current density, 
is almost purely azimuthal. 

The current density j follows the lines of the electric field E, and therefore behaves as 
sketched below: 



j, = o 




j ± =o- 

The total current can be evaluated by integrating the current density at r = a in region 



I: 



I = 2a r j r d9 = 2ao x f° E r {r = a)d9 « -2ka x f° d9 = -2ka x a = -2./-^^. (24) 
Jo Jo Jo y I — a 



In the present problem the total potential difference A0 is -2, so the effective conductivity 
is 

a= ^i = i— a (25) 

For a square checkerboard, a = 7r/4, and the effective conductivity is cr = *Jo\a2. It 
turns out that this result is independent of the ratio e^/cri, and holds not only for the corner 
region studied here but for the entire checkerboard array [S]. 
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